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Abstract. The two-dimensional a!- wave superconducting state of the high 
temperature superconductors has a number of different elementary exci- 
tations: the spin-singlet Cooper pairs, the spin 5 = 1/2 fermionic quasi- 
particles, and a bosonic 5 = 1 resonant collective mode, 4> a , at the an- 
tiferromagnetic wavevector. Although the <f> a quanta are strongly coupled 
to the gapped quasiparticles near the (n,0), (0,tt) wavevectors (the "hot 
spots"), they are essentially decoupled from the low energy quasiparticles 
near the nodes of the superconducting gap. Consequently, distinct and in- 
dependent low energy quantum field theories can be constructed for the 
(p a and nodal quasiparticle excitations. We review recent work introduc- 
ing a 2+1 dimensional boundary conformal field theory for the damping 
of the 4> a excitations by non-magnetic impurities, which is built on the 
proximity to a magnetic ordering transition at which the (p a condense; the 
results are compared with neutron scattering experiments. Photoemission 
and THz conductivity measurements indicate that the nodal quasiparticles 
undergo strong inelastic scattering at low temperatures; we propose that 
this is due to fluctuations near a quantum phase transition, and critically 
analyze candidate order parameters and field theories. 
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1. Introduction 

The description of high temperature superconductivity in the cuprate com- 
pounds has been a central problem at the frontier of quantum many body 
theory in the last decade. Although many anomalous properties have been 
observed in the normal state, both in the over-doped and under-doped re- 
gions, no theoretical consensus has emerged on their origin. Part of the 
difficulty is that there appear to be many competing instabilities and ex- 
citations as one cools down from high temperatures (T), and they are all 
strongly coupled to each other at intermediate T. 

However, simplifications do occur at temperatures T < T c , the critical 
temperature below which there is an onset of d-wave superconductivity. 
In this review we shall argue, on the basis of recent experimental obser- 
vations, that there is an important decoupling between different sectors 
of the excitation spectrum which carry a non-zero spin, and that this de- 
coupling allows development of tractable quantum field theories of the low 
energy excitations |3], |j . We will make quantitative predictions for the 
impurity-induced and intrinsic damping of these excitations and compare 
them to experimental results. 

Let us list the elementary excitations of the d-wave superconductor and 
nearby phases: 

(A) Cooper Pairs: The superconductivity is of course a consequence of 
the condensation of spin 5 = 0, charge 2e Cooper pairs. Below T c , the exci- 
tations of the phase of the condensate are responsible for the superflow, and 
for the plasmon excitations. In this paper, we will be primarily concerned 
with the damping of excitations which carry spin, and these couple only 
weakly to the phase excitations in a well-formed superconductor at low T: 
so we will neglect the phase excitations in the body of the paper. These 
phase excitations become more important near a T = superconducting- 
insulator transition, but we will not consider such a situation here. Above 
T c , phase fluctuations |^, |(| [7|, ||, |f| are surely important for the transport 
properties, and they also couple strongly to some of the fermionic quasi- 
particle excitations: we will briefly discuss this phenomenon further below. 

(B) S = 1/2 fermionic quasiparticles: These are the familiar Bogoliubov 
quasiparticles in a BCS theory of the superconducting state. Because of 
the d-wave symmetry of the order parameter, their energies vanishes at 
four nodal points in the Brillouin zone - (±if, ±K), with K = 0.3917T for 
optimally doped Bi2Sr2CaCu20g+,5 jlO|| . We will denote the fermionic ex- 
citations in the vicinity of these points by the Nambu spinors (see 
Fig [l] and further details in Section |3|). It is also interesting to consider 
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Figure 1. Brillouin zone of the high temperature superconductors at optimal doping 
(see e.g. Ref. @). The dashed line is the location of the incipient Fermi surface at 
intermediate temperatures: the ground state is a d-wave superconductor, and not a Fermi 
liquid, and so there is no sharply defined Fermi surface as T — > — the line is merely the 
location of smooth crossover in the momentum distribution function. The fermionic, 
S — 1/2, quasiparticles ^1,2 lie near the nodal points (±K, ±K) (with K w 0.39?r (!(]]) 
at which their excitation energy vanishes. The VP/, quasiparticles require an energy « A 
for their excitation and lie in the vicinity of the "Fermi surface" points (±0.18tt, ±7r) 
and (±7r, ±0.187r). The double-headed arrow at wavevector Q = (tt, tt) represents the 
bosonic, S — 1, resonant collective mode <f> a scattering fermions between two points in 
the Brillouin zone. Notice that the ^1,2 fermions are decoupled from the <f> a quanta. In 
contrast, the fermions couple strongly to the cj> a , especially in the vicinity of the "hot 
spots" denoted by the small filled circles. 



4 



the fermionic excitations near the (ir, 0) and (0, tt) points: here the pair- 
ing amplitude has its largest value and so there is a large energy gap, 
A, towards exciting the quasiparticles. We will denote these high energy 
quasiparticles generically by ^>h (see Fig[l|). Because of their large pairing 
amplitude, the ^>h quasiparticles couple efficiently to the phase fluctuations 
discussed above in (A), and are expected to have a rapidly decreasing life- 
time once the phase fluctuations proliferate above T c |]12[ . In contrast, the 
nodal quasiparticles, ^1,2, are in a region of vanishing pairing, and are es- 
sentially decoupled from the phase fluctuations: as we will discuss below, 
other mechanisms will be required to damp the ^12 quasiparticles. 



(C) S = 1 resonant collective mode: Neutron scattering experiments ob- 
serve a sharp resonance peak at an energy A res in the scattering cross 
section at the antiferromagnetic wavevector, Q [l4[ |l5|, [l6|, 17]. We 
will view this bosonic 5 = 1 resonant collective mode, <p a (a = x,y,z are 
the spin components), as that expected in a paramagnetic phase across a 
magnetic disordering quantum phase transition fTsif . (For the case where 
Q = (vr, 7r), the <j) a are real, while for incommensurate Q, the (ft a become 
complex; we will explicitly treat the commensurate case here, although the 
generalization to the incommensurate case is straightforward Q, and does 
not modify any of the scaling arguments, including the central result (^).) 
Our identification of <j) a is similar to the view that it is a S = 1 particle- 
hole bound state in a d-wave superconductor |l^, |2^, [2l|| . However, as we 
shall discuss in more detail below, appealing to the proximity of a quantum 
phase transition allows a systematic treatment of the strongly relevant self- 
interactions of the 4> a field. Such a collective mode has also been discussed 
in models with a special SO (5) symmetry [22], but we shall not appeal to 
symmetries beyond the usual SU(2) spin symmetry in our treatment. The 
coupling of the cp a to the fermionic quasiparticles in (B) is illustrated in 
Fig [l]. Momentum conservation allows a strong coupling between the 



and the 4>a in the vicinity of the so-called "hot spots" [11, 24, 25]. This 



coupling leads to strong mutual damping of ^ and <ft a above T c . Below T c , 
the same coupling is surely an important ingredient in the pairing of the 
fermions [26|: this conclusion is supported by the "shake-off" satellite 
peaks, separated by the A rcs ~ 40 meV, observed in the photoemission 
and optical spectra [^4], ^7|, A key ingredient in our discussion is also 
clear from Fig ||: the coupling between the low energy V&i a fermions and 
the 4>a mode is strongly suppressed by momentum conservation — moving 
a distance Q from a nodal point places one in a section of the Brillouin 
zone (e.g. near the point P in Fig ||) where the fermionic excitations cost 
over 100 meV. Remnants of the shake-off peaks just noted are also seen 
along the (1, ±1) directions p^] not too far from the point P, but this does 
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not qualitatively modify the low energy excitations within the shaded 
circles. 

We have now collected all the ingredients necessary to motivate our 
recent computations. 

In Section [2] we will discuss the T = broadening of the <p a collective 
mode by the substitution of a small concentration of non-magnetic impu- 
rities like Zn or Li on the Cu sites, and compare our theoretical results 
J||, |3| to experimental observations . It should be clear from the discus- 
sion above that, in the absence of such extrinsic broadening, it is possible 
for the <f) a resonance to be infinitely sharp at T = in a d-wave super- 
conductor. The (p a quanta couple strongly to the fermions, but these 
induce no damping as long as A res < 2A (A ~ 40 meV near optimal doping 



1 31 ] ) . Along the diagonals of the Brillouin zone, momentum conservation 
prohibits coupling to the gapless nodal fermions $12, and only allows cou- 
plings to quasiparticle excitations whose excitation energy exceeds A res . 

Section ^ will consider intrinsic T-dependent damping of the nodal 
fermionic quasiparticles by inelastic scattering [|], ||]. We have already dis- 
cussed the damping of the gapped ^ above. Below T c , photoemission ex- 



periments [32, p3, pTl[3j| observe negligible damping of the ^ fermions, 



and this is consistent with our considerations: the phase fluctuations in (A) 
are suppressed, while the 4> a mode in (C) leads to coherent pairing of ^>h 
fermions. Above T c , phase fluctuations proliferate, and their strong cou- 
pling to the *$>h is expected to lead to significant inelastic scattering of the 
additional scattering is also expected from the "hot-spot" coupling 
to the 4> a , and these expectations are consistent with experimental observa- 
tions. However neither of the fluctuations in (A) or (C) couple to the nodal 
fermions ^1,2^ so, with our present considerations we would conclude that 
the nodal fermions should be very sharp both below and above T c . The ac- 



tual experimental situation is dramatically different [10, 36] — these nodal 
quasiparticles have a large inverse lifetime, which decreases roughly linearly 
with T, and an imaginary component of a self energy which is roughly lin- 
early proportional to frequency, oj, for hu > k^T . Moreover, these damping 
rates change smoothly through T c , with little sign of the superconducting 
transition. We have to appeal to other inelastic damping mechanisms to 
explain these observations, and we will present a critical classification of 
candidates in Section ||[ 

2. Impurities and the S = 1 resonant collective mode 

Before discussing the effect of impurities, we state our model for the spin 
collective mode, 4> a , in the clean d-wave superconductor. A popular ap- 



proach in recent work [TL9L 20, Ell has been to compute the dynamic spin 
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Figure 2. (a) Particle-hole propagator which has a S = 1 bound state at the antiferro- 
magnetic wavevector representing the cj> a quanta, (b) and (c) Three- and four-point self 
interactions of the spin excitations. 

susceptibility of the underlying electrons and to identify this mode as a 
bound state pole near the antiferromagnetic wavevector: this is schemati- 
cally indicated in Fig ||a. However all such studies have so far neglected 
three- and four-point (and higher multi-point) self interactions of the spin 
excitations, which are schematically indicated in Figs ||b and [2|c respec- 
tively. One of the key points of our work is that it is essential to include the 
interactions in Figs ||b and |2|c in the low energy theory: these are strongly 
relevant perturbations, and in a sense, their effective strength is infinite 
(higher multi-point interactions can, however, be neglected). The simple 
results for impurity-induced damping we shall quote below rely on hyper- 
scaling properties, and these are a direct consequence of the self interactions 
in Fig |. 

So how does one obtain a tractable theory which includes the interac- 
tions in Fig ^ ? Our strategy is to appeal to the proximity of a T = 
magnetic ordering transition at which the excitation energy of cf> a , A res , 
vanishes. This is a quantum phase transition driven by the condensation 
of 4> a to a state with coexisting superconductivity and collinear spin den- 
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sity wave order. We will develop a theory for the quantum-critical point 
of this transition, and then use powerful field-theoretic methods to expand 
back into the region where A res is non-zero. The values of the interactions 
in Fig ^ will be universally determined by the underlying structure of the 
expansion in relevant perturbations of the critical field theory. Such an 
approach effectively reduces to an expansion in A rcs /J, where J is a mi- 
croscopic exchange constant, and the smallness of this ratio is our primary 
assumption. 

We shall use a theoretical model (which is supported by the 'pseudo-gap' 
phenomenology) in which the magnetic ordering transition occurs while the 
excitation energy of the ^ti remains non-zero, and so these fermions can be 
neglected in the critical theory of the transition. As discussed above, the 
^l s 2 fermions have little coupling to the magnetic excitations, and so an 
action for the transition can be expressed in terms of the cj) a alone. By anal- 
ogy with theories developed for the magnetic transition in insulators, and 
using general symmetry arguments [37], we can write down the following 



effective action for the bulk 4> a fluctuations in the d-wave superconductor: 



d 2 x / dr 



(i) 



the fermionic excitations with spin, ^1,2, have been integrated out, and 
for reasons already discussed, serve only to renormalize the values of the 
couplings in ([!]). The parameter c is the velocity of spin- waves in the or- 
dered phase, and s tunes the system between the two phases which lie on 
either side of a critical value s = s c . The quartic non-linearity, go, cor- 
responds to the interaction in Fig |2|c. How about the cubic coupling in 
Fig ||b ? This is implicitly accounted for in (|l|) in a manner we now de- 
scribe. By momentum conservation, if two of the particle-hole propagators 
in Fig ||b carry momentum Q, the third must carry momentum 2Q ~ 0, 
i.e., it represents the ferromagnetic spin component L a . This has allowed 



three-point couplings with the (fi a , including the kinematic term [37| in the 
action ~ ie a p 1 L a <ppd T 4>. T The L Q fluctuations are not critical, and after in- 
tegrating them out, one obtains renormalizations of terms already present 
in© & 

The magnetic properties of the bulk quantum phase transition described 
by of Sb have been worked out in some detail fl38| , |l8|| , and many aspects 
are in agreement with trends in NMR and neutron scattering experiments 
on the high temperature superconductors [39, fttj. Here, we will only need 



a few well-known scaling properties of the critical point of (|l|). Upon inter- 
preting third spatial dimension, can also represent the partition 
function of a classical Heisenberg ferromagnet in dimension D = 3 at finite 
temperature, and its Curie transition corresponds to the quantum-critical 
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point we are interested in. Both occur at a critical value r = r c , where there 
the functional integral over (ft a is invariant under the scale transformation 

x — > x/b 

T — > T jh 

<Pl b 3 - 1 ^^, (2) 

where b is a rescaling factor, and uh and tjh are known critical exponents 
of the D = 3 classical Heisenberg model. The last transformation in (||) 
represents the mapping of the composite operator <j)\, and its scaling di- 
mension is not simply twice that of cj) a because of corrections due to the go 
interaction in (|l]). 

We now turn to the effect of a dilute concentration of impurities. We 
will outline the central ingredients leading to our main result, and refer the 
reader to Ref. || for further details. Consider a single impurity at x = 0; 
by "impurity" we mean an arbitrary localized deformation in the vicinity 
of x = 0. One consequence of any such deformation will be a change in the 
value of s near x = 0, and this will lead to the following term in the action 

C J drtf>l(x = 0,r). (3) 

Under the scale transformation , we see immediately that £ has scaling 
dimension 1/z/f/ — 2 ~ —0.57, and is therefore irrelevant at the critical point 
of (|l]), and it has only weak effects on the bulk properties. To obtain a local, 
relevant, perturbation on the bulk fluctuations, we need to consider quan- 
tum mechanical effects associated with Berry phases. The Berry phases 
accumulated by the precession of spins in the host antiferromagnet cancel 



almost completely upon an average over the lattice sites [37]; however, in 
the presence of impurities it is entirely possible that this cancellation is 
disrupted, and a residual Berry phase of spin S (S must be an integer or 



half-odd- integer) survives |4l| , 42]. To account for this Berry phase we in- 
troduce a single unit vector n a (r) (n^,(r) = 1) representing the orientation 
of the net uncompensated spin, and the action 

S imp = isJdTA a (n)^f±, (4) 

where A a is a function of n a defined by e a ^dA^jdn^ = n a . This Berry 
phase is intimately connected to the fact that an external magnetic field 
will lead to a Curie susceptibility = S(S + 1)/3T from the impurity spin 
in the non- magnetic phase (in the absence of Kondo screening-see below) ; 
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this response is divergent as T — > and is a reflection of a (25 + 1)- 
fold degenerate level near the impurity. For the case of a non-magnetic Zn 
or Li ion replacing a magnetic 5 = 1/2 Cu ion in the high temperature 
superconductors, the above arguments on Berry phases strongly suggest 
that each such impurity should contribute a term like (0) with 5 = 1/2. 
This conclusion is supported by NMR experiments ^5|, and we will 

assume its validity in our discussion below. Other authors have modeled 



the Zn ion solely as a non-magnetic scatterer in the unitarity limit [46 



47, 48, |49|, |50[|. In such a model, the fermionic quasiparticles form quasi- 
bound states at the impurity sites at the Fermi level; we believe that after 
accounting for the strong local Coulomb repulsion at the impurity site, each 
bound state will capture only a single electron, and the low energy physics 
will then be described by (f§) with 5 = 1/2 (see also Ref @). 

We now need to couple the impurity degree of freedom, n a , to those of 
the host. The most important coupling is the simple linear term 

Sc = 7 J dTn a (T)<j) a (x = 0, t). (5) 

To compute the scaling dimension of 7 at the fixed point where the impurity 
and bulk degrees of freedom are decoupled, we note that if n a — > n a under 
the transformation (|2|), the impurity action (||) remains invariant. Under 
such a mapping, 7 has dimension (1 — t]h)/2 ~ 0.48. Unlike £, the coupling 
7 is therefore relevant at the 7 = fixed point, and plays a central role in 
the main results presented below. We can also imagine a Kondo coupling, 
Jk between the spin, n a (r), and the host fermions ^1,2, ^h- However, the 
fermionic, single particle density of states vanishes at the Fermi level, and 
this dramatically reduces the possibility of Kondo screening of the impurity 
spin: with particle-hole symmetry, there is no Kondo screening even upto 
Jk = 00, while without particle- hole symmetry, the spin is screened only 



above an appreciable threshold value of Jk |>l|, |52|, |53| |5j. We will assume 
that no Kondo screening has occurred over the experimentally relevant 
temperature range. 

The above arguments suggest that a theory of the impurity spin dy- 
namics will emerge from a complete renormalization group analysis of 
Sb + 5i mp + S c . This has been carried out in Ref ||, and the final results 
are quite simple: the bulk phase transition at s = s c is the only critical 
point, and (s — s c ) remains the only relevant perturbation at this critical 
point; both g and 7 approach fixed point values g* and 7* (related phenom- 



ena were noted earlier in simpler models 56, |57||). There is no separate 
critical point associated with the impurity degrees of freedom, as is often 
the case in the theory critical phenomena on boundaries |p8| . A remarkable 
consequence of this is that the single energy scale, A res , which character- 
ized the dynamics of the paramagnet in the host system (37]], is also all 
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that is needed to completely characterize the dynamics in the vicinity of 
the impurity. Of course, the quantized number S in (Q) also influences the 
values of the various universal scaling functions. 

Now consider a dilute concentration of impurities, ni mp , each described 
by the analog of (|j) and (|5|), placed at random locations in the cf-wave 
superconductor. We will answer the following key question, relevant to the 
neutron scattering experiments on Zn doped YBa2Cus07 by Fong et al. 



1 3C ] . At what energy scale, T, does the S = 1 resonant pole at energy A res 
get broadened by the impurities ? The implication of the renormalization 
group arguments above is that, for small A res /J, T is universally determined 
by the only dimensionful parameters available to us: ni mp , A res , and the 
velocity c. Making the mild assumption that, for small ni mp , T must be 
linearly proportional to ni mp (this is supported by explicit computations 
Q), simple dimensional analysis of the length and time scales allows us to 
conclude 

r = C s ^n imp , (6) 

^rcs 

where C$ is a universal number. All corrections to (||) will be suppressed by 
positive powers of A res /J. Notice also the inverse dependence on the small 
energy scale A res : this is an indication of the strong effect of the relevant 
coupling in (||). For an impurity with S = 0, the simplest allowed coupling 
is £, and its irrelevance implies Cq = 0, and the corrections just noted will 
be the leading contributions. We have estimated C1/2 in a self-consistent 
non-crossing approximation and obtained Cu2 ~ 1. 

It is useful to rewrite our main result @ in a different manner: 

p 

-r n imp £ 2 , (7) 

'-ires 

where £ = hc/A Tes is a correlation length. So if we imagine a "swiss cheese" 
model |59| where each impurity makes a hole of radius £, the inverse Q of 
the resonance is of order the fractional volume of holes in the swiss cheese. 

The numerical predictions of (||) are in good agreement with the obser- 
vations of Fong et al. []30"| . We use ni mp = 0.005, A res = 40 meV, and the 
spin-wave velocity in the insulator he = 0.2 eV, and obtain T = 5 meV. 
This compares well with the observed value of 4.25 meV. We have also pre- 
dicted detailed lineshapes for the impurity-induced broadening, and these 
will hopefully be tested in future, higher precision experiments. 

3. Inelastic damping of the nodal quasiparticles 

As we noted at the end of Section |], recent experimental observations 
|10|. |36| of a short lifetime ~ h/k^T, and a large frequency-dependent self 
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energy, for the nodal quasiparticles both above and below T c , are puzzling 
in the light of the very weak coupling between the \&i 5 2 and both the <p a and 
phase fluctuations. Consequently, we have to appeal to a separate decou- 
pled sector of low energy fluctuations to explain this anomalous damping 
1,1- 

A natural way of obtaining inverse lifetimes of order ksT/h, and self 
energies of order u>, is to assume that the system is in the quantum-critical 
region of a T = quantum phase transition for which the ^12 fermions are 
central critical degrees of freedom [41, (see Fig|||); the quantum-critical 



point should be described by an interacting quantum field theory below its 
upper critical dimension, so that universal low-energy fluctuations dominate 
the interactions. A further constraint is that the critical fluctuations should 
be decoupled from the fermions, as these remain undamped below T c : 
so the new low energy mode associated with the onset of state X in Fig |3] 
couples strongly to the but not the $>h- It will turn out that these 
constraints are rather difficult to satisfy, and lead to an essentially unique 
identification of the state X in Fig ||. 

Clearly, the magnetic transition in Section [2] cannot be the required 
transition because the ^^2 fermions are innocuous spectators of its critical 
field theory. We will list below many of the order parameters that have 
been considered in the literature in the last decade, and discuss whether 
they satisfy the requirements we have imposed on the state X. 

Before we embark on this, let us recall the effective action for the nodal 
quasiparticles in the <i-wave superconductor well away from the quantum- 
critical point, r S> r c . We denote the components of the electron an- 
nihilation operator, c a , in the vicinity of the four nodal points (K,K), 
(-K,K), (-K,-K), (K,-K) by f la , f 2a , f 3a , ha respectively, where 
a = j, I is the electron spin component. The 4-component Nambu spinors 

are *i = (fi a ,£abfl b ) and ^ 2 = (f2a,£abfl b ) where e ab is an antisymmetric 
tensor with e-m = 1. The action is then 

f d d k 

Sy = / — TjT^f} (-iu) n + v F k x T z + v A k y T x )V 1 

+ J S^T^l (-^n + V F k y T z + V A k x T X ) tf 2 . (8) 

Here r a are Pauli matrices which act in the fermionic particle-hole space, 
k X)V measure the wavevector from the nodal points and have been rotated by 
45 degrees from the axes of the square lattice, and vf, v a are velocities. The 
action 5<i> has a scale invariance which will be useful in our considerations 
below 

x — ► x/b 



12 



T 



o 



T 



7~, 



QUANTUM 
N n CRITICAL / 



x 



\ 



/ 



\ / 
\ / 

superconducting X\ / d-wave 



state X 



superconductor 



Figure 3. Finite temperature (T) phase diagram in the vicinity of a second order 
quantum phase transition from a d-wave superconductor as a function of some parame- 
ter in the Hamiltonian, r (which is possibly, but not necessarily, the hole concentration 
5). Superconductivity is present at temperatures below T c , and the superfluid density is 
non-zero on both sides of r c . The state X is characterized by some other order parameter 
(in addition to superconductivity) which vanishes above a temperature Tx- We discuss 
a number of possibilities for the state X in the text-only one naturally satisfies the re- 
quirement of leading to an inverse lifetime ~ ksT/h for the nodal quasiparticles \ti,2 
in the quantum-critical region, and negligible damping of the quasiparticles below 
T c : the (d x 2_ y 2 + id xy )-wicve superconductor. Our computations for the quantum-critical 
point are carried out below T c , but the results should also apply above T c as long as 
the quantum-critical length ~ j 1-1 / 3 remains shorter than the phase coherence length. 
We emphasize that we are not requiring the high temperature superconductors to have 
(d x 2_ y 2 + id X y) order in the ground state (although it is permitted): the coupling r 
could be larger than r c , but it should be close enough that the system enters the quan- 
tum-critical region at some low T. 



fi, 2 (x,r) 



r/b 

bV 1>2 (x/b,T/b), 



(9) 



where 6 is a rescaling factor. We can illustrate the power of scaling ar- 
guments by showing how this transformation allows us to quickly deduce 
damping produced by ordinary screened Coulomb interactions (the screen- 
ing is performed by the Cooper pairs in the condensate). These interactions 
lead to couplings like 



d 2 x 



(10) 
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and it is easy to see that under (||), v has scaling dimension -1. Therefore v 
is irrelevant, and perturbation theory in v should be reliable. The fermion 
self energy, £j, acquires an imaginary part at order v 2 ; using the fact that 
under (^) the frequency and momentum dependent self energy transforms 
as X/ — > the result (|To| ) immediately leads to ImSj ~ v 2 T 3 (and 

similarly for the u dependence). This damping rate appears far too weak 
to explain the experimental observations. 

We now turn to considerations of quantum critical points, and list vari- 
ous plausible candidates for the state X in Fig |J| We will only consider order 
parameters for X which break simple underlying symmetries of the Hamil- 
tonian: the symmetry of the square lattice space group, time-reversal, and 
spin rotation (the last we have already discussed above). More complicated 
transitions with non-local order parameters and deconfinement transitions 
are also possible, but we will not consider them here J6C|]. 



Staggered flux phase: Many investigators 32, 63, 64, 65, 56, 67] have 



considered the possibility of d-wave superconductivity coexisting with a 
staggered distribution of orbital currents (or an "orbital antiferromagnet" ) . 



This state X is characterized by the expectation value [66| 



( c k+Q,a C k,a) = i(t*(cOSk x - COS ky) , 



(11) 



where, as before, Q = (vr,7r), and is a real order parameter. Momentarily 
neglecting the fermionic excitations, we can, just as for ([lj), write down an 
effective action for eft fluctuations purely on symmetry grounds: 



d 2 x 



dr 



-((5 r 0) 2 + c 2 (V :I 



? + r6 2 



(12) 



(without strict particle-hole symmetry, a first order time derivative term 
is potentially allowed, but the only possible relevant term, (pd T (j), is a total 
derivative (l/2)d T (p 2 ). As is well known, <S^ describes the phase transition 
in the Ising model in D = 3 spacetime dimensions. At the critical point 
r = r c , Sb is also invariant under the extension of (||) to the analog of the 
scale transformations in (||) 

^ ^ 63-1/^2 f (13) 

where now vj and rjj are the exponents of the D = 3 Ising model. To 
decide if this purely Ising description of the transition is correct, we have 
to test its stability to a coupling between the <p and the ^1,2- From (pd|) 
we see that cj) carries momentum Q, and so the momentum conservation 
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constraints upon its coupling to the fermions are identical to those of the 
spin mode c/) a in Fig |]: unless K = tt/2, there is no linear coupling to 
the fermionic excitations which is linear in (f>. For general K, the simplest 
allowed coupling is 

w J d 2 x J dr^ 2 *|r 2 *i, (14) 

and similarly for ^2- Using (13), we deduce that the scaling dimensions 
of w is l/uj — 2 ~ —0.41: consequently w is irrelevant, and the critical theory 
for the transition is ( |l^ ) alone. The fermions V&i^ are not part of the critical 
theory, and their inverse lifetimes can be estimated by perturbation theory 
in w. Using the scaling dimension of w above, we deduce that the fermionic 
self energy has the following T dependence in the quantum-critical region 

Of Fig H ImS/ ~ w 2 T 5-2/„j _ W 2 T 1.83 ( and ImS/ ^ ^2^1-83 for ^ > 

fcflT). This is a super-linear power, which does not appear compatible with 
experimental observations. Finally, we note that the special case K = n/2 
has also been considered in Ref 0: then a coupling term linear in <j> is also 



allowed [66|, but its contribution to £j vanishes with an even higher power 
of T. 



Charge Stripes: We consider the onset of a charge density wave in a cf-wave 
superconductor, a transition to a state X defined by the order parameters 

( C k+G x ,a C k,a) = $x ] (4+G„a C M) = $y (15) 

where G x = (G,0), G y = (0, G) are the ordering wavevectors, and &x,y 
are complex order parameters. Again, constraints from momentum conser- 
vation are rather severe. Unless G = 2K, there is no coupling between the 
Q x ,y and the fermions. Existing experimental observations of charge 
stripe formation easily satisfy G ^ 2K. Under these conditions, the damp- 
ing of the ^1,2 from the critical charge fluctuations can be estimated as in 
the staggered-flux case above: the simplest allowed couplings, as in (PH), 
are ~ |$ x | 2, i4t^i etc., and we obtain ImSj ~ (max(w, T)) 5 ~ 2 /^. Any 
reasonable model of the critical theory of the & x ,y 01 has v > uj, and so 
the fermion damping is rather weak. The special case G = 2K has also 
been analyzed in Refs ffl: it does yield fermion damping compatible 



with experimental observations [10, 36|, but, as we have already noted, this 



mode-locking of the charge stripe and fermionic nodal wavevectors is not 
supported by experiments. 



d + is superconductivity: Next consider a time-reversal symmetry break- 
ing transition in which the Cooper pair wavefunction in state X acquires a 
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small s-wave component, but with a relative phase factor ±7r/2 [68]: 



(cfcjC-fcj) = Ao(cos k x — cos k y ) + i<f>(cos k x + cos k y ). (16) 

The order parameter is again a single real field 0. On general symmetry 
grounds, we expect the effective action of the <f> fluctuations to also have the 
form (^). However, now there is an efficient coupling of cf> fluctuations to 
the nodal fermions, which is not preempted by momentum conservation. It 
is evident from (^) that 4> fluctuations Andreev scatter fermions with mo- 
menta k and —k. This scattering can occur between the nodal points (K, K) 
and (—K,—K) of the fermions, and similarly for \p2> it is represented 
by the allowed coupling: 



J d 2 xdr [a o (V^l-i + ^r^ 2 ) ] • (17) 



Now, computing the scaling dimension of the coupling Ao under the trans- 
formations (|9|), (|l~3|), we observe a crucial difference from the two cases 
considered so far: the coupling Ao has dimension (1 — rjj)/2 « 0.48, and 
is therefore relevant, and the Ao = fixed point is unstable. So the crit- 
ical theory strongly couples the <\> and fluctuations, and a complete 
understanding requires a more detailed renormalization group analysis of 
S\jj +<S0+«S$0. This has been discussed in Ref. [§]] (and for a similar model in 
a different physical context in Ref. (69]): we will not discuss this here apart 
from noting that both the non-linearities, g and Ao approach fixed-point 
values, and there is only one relevant perturbation, (r — r c ), at the inter- 
acting critical point. Under these conditions, strong scaling applies, and 
the lifetimes of excitations of the <j> and quanta are of order H/kBT in 



the quantum-critical region of Fig y [41, 37]. So for the case in which X 



is a d + is superconductor, the relaxation of the nodal quasiparticles V&i^ 
appears to be in good accord with experimental observations. However, 
one significant discrepancy remains: the s-wave order parameter couples 
strongly to fermions in all directions, and so will also couple to the 
quasiparticles (see Fig |lj) . The gapped ^ quasiparticles will easily radiate 
many of the low-energy (ft quanta, and acquire an appreciable width: this 
is not in accord with photo-emission experiments in which, as we noted 
earlier, the quasiparticles become sharp below T c . 

d x 2_ v 2 + id xy superconductivity: Finally, we consider another case in which 



X breaks time-reversal, but now by acquiring a small d xy component [70|. 
We will see that this case is very similar to the d + is case discussed above, 
but it also succeeds in very simply and naturally resolving the discrep- 
ancy with width we have just mentioned. The order parameter for this 
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J d 2 xdr [a o (*\t v Vi - ^ 2 T y ^ 2 ) } • (19) 



transition, replacing fllq), is 

(cfc|C_fc|) = Ao(cos/c x — cos k y ) + icp sin k x sin k y . (18) 

The <^ fluctuations again Andreev scatter fermions between k and —k, and 
their coupling to the Vl/i 2 fields has the form (replacing ([FT 

Notice that the only difference from ( |17|) is the relative sign of the \Pi and ^2 
terms: this is because the coefficient of <fr in (O) changes sign between the 
two pairs of nodal points, unlike the case in tfPjj). The consequences of (^) 
are essentially identical to those of (0): the coupling Ao approaches a fixed 
point value, and this leads immediately to a lifetime ~ fi/ksT for the nodal 
quasiparticles in the quantum-critical region. Moreover, the fermions do 
not couple to the </> fluctuations: the coefficient of <j> in (^) vanishes along 
the line between (ir, tt) and (tt, 0) (and also between (0, 0) and (jr, 0) and 
other symmetry-related lines), and so the fluctuating d xy component of 
the pair wavefunction does not lead to appreciable broadening of the *$>h 
quasiparticles. Therefore, if the state X is a d x i_ y i + id xy superconductor, 
the dynamics of the quantum-critical region are very naturally in accord 
with the constraints described at the beginning of Section ||. 

Quite apart from the motivation provided by the above analysis, there 
are some appealing independent reasons [ff(| for suspecting that the d-wave 
superconductor may be on the verge of an instability to a d x 2_ y 2+id xy state. 
The structure of the incipient Fermi surface in Fig [I] indicates that there is 
significant second-neighbor hopping on the square lattice. Accompanying 
this there should be a corresponding second- neighbor exchange, J2. Just 
as the first neighbor exchange, J±, prefers d x 2_ y 2 pairing, the J2 exchange 
will induce d xy pairing; the realistic case with both J\ and J% non-zero 
should therefore prefer the intermediate d x 2_ y 2 + id xy state, with relative 
phase of ±7r/2 ensuring that a gap opens over the entire fermion spectrum. 
Alternatively stated, the ^ fermions are already strongly paired in the 
d x 2_ y 2 state, while the fermions are essentially unpaired; the system 
will try to lower its energy by pairing the fermions, and this is most 
efficiently done by an additional d xy component to the pair wave-function. 
An additional s component would also do the job, but has the disadvantage 
of also deforming the already optimal pairing of the and so is not as 
efficient. Our mean-field calculations Q on models with J\, J2 both non- 
zero are consistent with these expectations. We can therefore identify the 
coupling r in Fig ^ as r ~ J\ / J 2 ■ 

It is also important to note that the instability from a d x 2_ y 2 super- 
conductor to a d x 2_ y 2 + id xy superconductor occurs below a finite value of 
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the coupling r = r c . This is to be contrasted with pairing instabilities of a 
Fermi liquid, which occur at infinitesimal attraction, and so the analog of 
the effective action ( |l~2| ) has a logarithmic dependence upon the order pa- 
rameter. However, when the parent state is a d x 2_ y 2-wave superconductor, 
the vanishing density of states at the Fermi level removes the usual BCS 
log divergence, and a finite attraction is required for further pairing in the 
d xy channel. Only such a finite- coupling quantum phase transition can be 
described by an interacting quantum field theory with hyperscaling prop- 
erties, and which leads to a T > quantum-critical region with lifetimes of 
order h/ksT. 

Further tests of the above scenario will be provided by computations 
of transport properties, including the optical conductivities and the Hall 
coefficient, in the quantum-critical region of Fig || for the case where X is 
the d x 2_ y 2 + id xy superconductor: these are currently in progress. 
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